Ratchet transport of interacting particles 
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We study analytically and numerically the ratchet transport of interacting particles induced by a 
monochromatic driving in asymmetric two-dimensional structures. The ratchet flow is preserved in 
the limit of strong interactions and can become even stronger compared to the non-interacting case. 
The developed kinetic theory gives a good description of these two limiting regimes. The numerical 
data show emergence of turbulence in the ratchet flow under certain conditions. 



PACS numbers: 05.60.-k, 47.61.-k, 72.40,+w 

I INTRODUCTION 

For systems without spatial inversion symmetry the 
appearance of directed flow of particles induced by a 
time-periodic parameter variation with a zero-mean force 
is now commonly known as the ratchet effect (see re- 
views [1, 2, 3]). This phenomenon is ubiquitous in na- 
ture so that such flows appear in a variety of systems 
including asymmetric crystals [4, 5] and semiconductor 
surfaces [6] under light radiation, vortexes in Josephson 
junction arrays [7], macroporous silicon membranes [8], 
microfluidic channels [9] and others. A significant in- 
crease of interest to ratchets is related to the experimen- 
tal progress in the investigation of molecular transport in 
biological systems like proteins characterized by asymme- 
try and non-equilibrium [1, 2, 3]. At the same time the 
nanotechnology development allowed to fabricate artifi- 
cial asymmetric nanostuctures with the two-dimensional 
electron gas (2DEG) where it has been shown that in- 
frared or microwave radiation creates a ratchet transport 
[10, 11, 12, 13]. The theoretical studies predicted that the 
directionality of ratchet flow in such systems can be con- 
trolled by the polarization of radiation [14, 15, 16, 17, 18] 
that has been confirmed by recent experiments with a 
semi-disk Galton board for 2DEG in AlGaAs/GaAs het- 
erojunctions [19]. 

Till present the theoretical studies of ratchet trans- 
port have been done mainly for non-interacting particles 
[1, 2, 3, 4, 5, 14, 18]. However, in many systems the inter- 
actions between particles are of primary importance like 
for example for microfluidic channels [9], 2DEG nanos- 
tructures with strong electron-electron (e-e) interactions 
at a large r s parameter [20], granulated materials [21] and 
one-dimensional Luttinger liquids [22]. On a first glance 
it seems that a strong scattering between particles should 
suppress the ratchet transport. But on the other hand 
the local conservation of momentum of particles indicates 
that even in presence of strong interactions the ratchet 
flow still should exist. The investigation of the proper- 
ties of ratchet transport for interacting particles in two 
dimensions is the main aim of this paper. The theory 




FIG. 1: (color online) Geometry of asymmetric scatterers ori- 
ented in (x, i/)-plane: cuts with elastic (left) and diffusive 
(right) sides; elastic semidisks; liner-polarized force F has an- 
gle 6 in respect to :r-axis. 

developed is based on the kinetic approach used in [18] 
extended to the case of strong interactions. The theory 
is compared with the extensive numerical simulations of 
ratchet transport of interacting particles in asymmetric 
structures. The model description is given in Section II; 
the analytical theory based on the kinetic equation is de- 
veloped in Section III; the numerical results are presented 
in Section IV and the discussion is given in Section V. 

II MODEL DESCRIPTION 

The interactions between particles are treated in the 
frame of the mesoscopic multi-particle collision model 
(MMPCM) proposed by Kapral (see e.g. [23]). This 
method exactly preserves the total momentum and en- 
ergy of particles colliding inside each of N ce i collision cells 
on which the whole coordinate space with N particles is 
divided. In this method the collisions inside cells are 
modeled by rotation of all particle velocities in the mov- 
ing center of mass frame of a given cell on a random angle 
after a time tk- To equilibrate the whole system of in- 
teracting particles in presence of external monochromatic 
driving force F cos Lot we use the Nose-Hoover thermostat 
[24] which drives the system to the Boltzmann equilib- 
rium with a temperature T — mvT 2 /2 on a relaxation 
time th- Such a combination of two methods for sys- 
tems with interactions and ac-driving has been already 
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FIG. 2: (color online) Polarization dependence of the average 
ratchet flow Vf in the flashing cuts model with (top panel) and 
without (bottom panel) interactions; diamonds and circles 
show numerical data for vf ]X and Vf lV components, curves give 
the fits of data (see text). The system parameters are: N = 
10 4 , N ce i = 100 x 100 inside the periodic space domain Rx R 
with v t th/R = 2.4, t c /t h = 0.45, t k /t h = 0.02, lot h = 3, 
FvTr c /T = 0.64 for the top panel and same parameters for 
the bottom panel but tk /th = oo and impurity scattering is 
added with Tim/rn = 0.5; total integration time is i/rn ~ 
10 3 . 



used in [25]. As in [18] the asymmetry appears due to 
asymmetric scatterers having form of vertical cuts with 
diffusive (right) and elastic (left) sides (cuts model) or of 
elastic semi-disks of radius Td (semi-disks model) placed 
in a periodic square lattice of size Rx R. The system ori- 
entation geometry and two types of scatterers are shown 
in Fig.l (see also [18] and Fig. 5 below). In the cuts 
model it is assumed that the scattering on cuts takes 
place instantaneously at random moments of time which 
have a Poisson distribution with time scale r c . This cor- 
responds to the case of flashing cuts model (instanta- 
neous appearance of cut at some moment of time) which 
is slightly different from the case of static cuts randomly 
distributed in the plane (both cases were discussed in 
[18]). As in [18], in absence of interactions an effective 
impurity scattering is added with the scattering time Tj m . 
The monochromatic force is polarized as it is shown in 
Fig. 1 with F = F(cos9, sin#). Here, we present numeri- 
cal results only for the semidisks model and the flashing 
cuts model, which is rather convenient for numerical sim- 
ulations, but in the analytical treatment we also consider 
the static cuts model. 

The results of numerical simulations for the polariza- 
tion dependence of the ratchet flow in the flashing cuts 
model are shown in Fig. 2. In absence of interactions 
the results are well described by the theory [18] with the 
fit dependence Vf/v? — 6(— cos(20), 2sin(20))/2 where 
b = 0.0064(Ft> T r c /T) 2 « 0.8b th and b th is the theory 
value (see Eqs.(9),(41) in [18]). For interacting particles 
the fit gives the dependence Vf/Vr = b; mt {— a\ cos 2 9 + 
a,2 sin 2 9, sin(20)) with bi nt /b — 2.7 and a\ = 0.10, 02 
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FIG. 3: (color online) Same as in Fig. 2 for the semi-disk 
model with R/rd = 4, Frj/T — 0.15, ljth = 1, effective 
Tc/th ~ R 2 /(2rdVTTH) = 0.85, other parameters are as in 
Fig. 2. 

0.29. In presence of interactions the ratchet flow appears 
even after polarization averaging. The results for the 
semi-disks model are shown in Fig. 3. Without interac- 
tions the data are satisfactory described by the theoret- 
ical dependence v f /w T = b(— cos(26>), sin(20)) with the 
fitting value b = 0. 24 (Fr d /T) 2 w 0Ab th and the theoret- 
ical value bth of [18] (see Eq.(42) and discussion there). 
In presence of interactions the polarization dependence 
of the flow is qualitatively changed: the component v y is 
enhanced by a factor 8 and v x remains negative for all 9 
showing signature of 4th 9— harmonic (Fig. 3, top panel, 
curves are drown to adapt an eye). 



Ill ANALYTICAL THEORY 

The numerical simulations are based on the dynamical 
description of motions of many interacting particles. To 
obtain an analytical description of the ratchet transport 
we use the kinetic equation approach valid for systems 
with developed chaos and rapid decay of correlations. 
The validity of the kinetic equation requires rare colli- 
sions with asymmetric scatterers (antidots) and random- 
ness of scattering events. Under such conditions the ki- 
netic equation can be applied for comparative study with 
the numerical data even if the numerical simulations are 
done for a deterministic system with a periodic lattice of 
semi-disks of relatively large size. 

The symmetry of the system determines the ratchet 
flow which mean velocity v/ is quadratic in the ampli- 
tude of the ac-force F(i) = Re(Fe~ iwt ). Therefore, the 
flow velocity can be described by the phenomenological 
expressions 

Vf.x — &xxx I F x I -(- Q?a;yy|-fy| : ^f,y — 2R,e(cVy ;C y F x Fy ) ■ 
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The tensor components a xxx , 



and He(a yxy ) de- 



termine the response produced by a linear-polarized 
monochromatic force (ImF = 0). In absence of inter- 
actions (see [18]), for the linear polarization along x or y 
axes the mean flow is directed along x-axis; the current 
in y direction appears for tilted linear-polarized force. 

We also note that for the elliptically-polarized force 
with ImF 7^ there exists also a circular ratchet effect 
determined by the product of lm(a yyx ) and lm(F x F*) 
but we will not consider this effect here. 



Kinetic equation 



The kinetic equation in the momentum space p reads 



(i) 



where in the case of microwave field, E(i) is the electric 
field E(£) interacting with electron gas F(t) — eE(t), e is 
the electron charge. The collision operator I = I e i + I ee 
contains the operator of elastic collisions (including im- 
purities and scatterers (or antidots)) I e i and interparticle 
(electron-electron or e-e) collisions I ee . 

The integral of elastic collisions with scatterers and 
static impurities reads as 

iei(u) = e^pp'/p' = Et H/ (p , 'P)/p' - ^(p.p')U 

p' p' 

(2) 

where Q pp > is the kernel of the operator I e i and W-^p', p) 
is the probability of the transition from p' to p. 

The interparticle collisions operator (e-e) operator is 

Iee{f) = ^2 H VtPLP'+P'i ( 3 ) 

Pl,p',p'l 

x 5(e p + e pi - e p , - e pi )u p - P ' 

x {/p'/ pi (l - / P )(l - / Pl ) - / P / Pl (l - / P ')(l - / Pi )} • 

Here S is the sample area, Uk is the Fourier transform of 
e-e-interactions. 

Interparticle collisions satisfy the conservation of the 
total momentum of gas. Due to the Galileo invariance the 
action of the collision integral on the equilibrium distri- 
bution function with shifted argument iee/ p + a vanishes 
for any a. Expanding by a we have: 

7 ee (/ p 0) ) = 0, 4(a5 p /(°)) = 0, 
^ e (a9 p /(°) * a9 p ,/(°)) + /LK«^, w / (0) ) = 0. (4) 

We use the following notations for the first and the second 
variations around equilibrium: 5I ee (f) = I' ee (Sf) (linear 
operator), S 2 I ee (f) — I' e ' e (Sf *Sf) (bi-linear operator, as- 
terisk denotes integration with two functions of different 
arguments). 



The ratchet flow is generated by the anisotropy of col- 
lisions. This anisotropy is constructed artificially due to 
asymmetric form of oriented scatterers. As theoretical 
models we considered cases of fixed oriented anisotropic 
scatterers, namely cuts and semidisks. The model of 
static cuts is analytically solvable [18] but has a disad- 
vantage since it leads to a divergence due to electrons 
moving along the mirrors. Even if this divergence can 
be regularized by an isotropic impurity scattering such a 
property is not very convenient. Due to that it is useful to 
use a modified model of flashing cuts which does not have 
such divergence. In this model at any moment a particle 
can meet a scatterer with a constant probability indepen- 
dent of its velocity and direction of motion; after collision 
the particle equi-probably scatters into any angle of the 
right semicircle if it collides from the right semicircle and 
is mirror-reflected if it collides from left semicircle (see 
Fig.l). Such a model of flashing cuts gives a significant 
simplification for analytical and numerical studies. 

The corresponding transition probability in these mod- 
els are given by (see also [18]): 

4tt 2 

W(p',p) = —w(<p',<p)6(e p -e p ,); (5) 



with 



((f',ip) = — COS if' 6>(cOS if')S((f' + tp — 7r) (6) 



— — cos (p' cos <p 9(cosLp)9(— cosy/) (static cuts), 



w(ip', if) = — cos f' 6>(cos Lp')S(ip' + tp — 7r)+ 

I|sin(^^)|[^-^(-V-^) (7) 
+ 6(<p' — <p)8(tp + <p')] (semi-disks), 



w((p', tp) — — [#(cos ip')5((p' + <p — 7r) + 

— Q (cos ip)6(— cos <p')\ (flashing cuts). (8) 

Here <p is the polar angle of electron momentum (— tt < 
p < 7r), r c (e) is the characteristic scattering time on 
asymmetric scatterers, 8(x) is the Heaviside function. 



Linear response 

We consider the limit of high rate of interparticle scat- 
tering exceeding the rate of elastic collisions. At the same 
time the interactions preserve the total momentum and 
in isotropic media do not affect the momentum relax- 
ation. This is not the case for an anisotropic medium 
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where the interactions indirectly lead to the momentum 
relaxation due to the conversion of the first angular har- 
monic of the distribution function f p to higher harmonics 
produced by the anisotropic scattering. In particular, it 
is generally excepted that in an isotropic medium with 
closed Fermi surface the e-e scattering does not affect 
the conductivity. Nevertheless, in the considered case 
of anisotropic medium e-e collisions indirectly affect the 
momentum relaxation rate. This action is realized due 
to the conversion of the first angular harmonics of the 
distribution function to higher harmonics produced by 
the anisotropic scattering. As a result, the conductiv- 
ity becomes temperature dependent in the temperature 
range when the e-e relaxation time is comparable with 
the elastic relaxation time. 

At first we consider the linear response to the electric 
field using the expansion / = /(°) + + f^> + ... in 
small driving force F. The linearized kinetic equation 
can be written in the form if {1) (t) = Re(/i, 1) e^*)): 



and 



iufV+F u d p fW=lV(fW), 



(9) 



where the collision operator contains the elastic collisions 
with anisotropic scatterers determined by I e i and inter- 
particle or e-e collisions determined by I' ee : 



(10) 



The formal solution of Eq. (9) in the first order of alter- 
nating force is 



/( 1 ) = (^ + 7«)- 1 (F^ p )/(°). 



In the case of weak e-e interaction I' 



be 



(11) 



clcd. 



In the opposite limit of strong e-e scattering the formal 
parameter describing I ee is large. Having in mind Eq. (4) 
we see that the inverse operator (ui + F 1 ^)^ 1 can be found 
by a projection on the subspace of the Hilbert space of 

the basis functions Vi = g p . /II g p . II corresponding to 

zero eigenvalue of the operator I' ee . Thus the operator 
I e i is replaced by its projection, while I' ee can be can- 
celed. The resulting tensor of conductivity of e-charged 
particles with density n e reads 



(Tij (uj) = 



e 2 n R 



m 1 — iujTi 



(12) 



where Ti are relaxation times of the first harmonics of the 
distribution function related with the projected operator 
of elastic collisions: 



— = -^^*(p)<9pp'^i(p')- 



(13) 
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8 + 7T 
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(for static cuts), 
(for semi-disks), 



a x = — H — -, a y = — (for flashing cuts). (14) 

The quantity f c is determined by gas statistics: 
1 _ ^de{f^f{e/r c {s)) 

where prime notes the derivative over the energy e. 

In the case of static cuts or semi-disks l/r c (e) oc e s 
with s = 1/2. So one can write f c = t c (sf) (strongly 
degenerate Fermi case) and f c = 4y / 2/7rT c (T)/3 (Boltz- 
mann case); s — 1/2 for fixed obstacles and s = for 
flashing cuts (in this case r c (e) = const). 

The physical origin of Eqs. (12) and (13) is a very 
quick relaxation of higher angular momenta harmonics 
as compared to the first harmonic relaxation. As a re- 
sult the conductivity has different values at low temper- 
ature, when T ee 3> T e i and at high temperature when 
T ee <C Tei- In both limits the conductivity does not de- 
pend on e-e interaction, but has different values. In the 
case of the Fermi distribution the conductivity changes 
from low temperature value where r ee ^> T e i to high tem- 
perature value where r ee <C r e ;. We should emphasize 
that the transition between these two values is ruled by 
the ratio T ee /r e / rather than by the ratio of temperature 
T to the Fermi energy Ep. The transition temperature 
To can be estimated by equating e-e relaxation time to 
the relaxation time given by elastic scattering. In clean 
samples with high mobility the transition corresponds to 
a rather low temperature To ~ Ep/a yj\p/l pi where 
a = (e*) 2 /hvp is the dimensionless e-e interaction con- 
stant, Xf and vf are the Fermi wavelength and veloc- 
ity l p is the elastic mean free path. For Ep = O.OleV, 
X F w lOnm, a = 0.5, l p - 10~ 4 cm, T - 10K. 

From Eq.(12) one can write the expression for ratio of 
static conductivities <J vy /a xx : 

0~yy/&xx = Ty/Tx = CLx/dy- (15) 

In case of flashing cuts this ratio is equal to 3 + 8 /ir 2 w 
3.81 (sec Eq. (14)). For such scatterers the problem of 
linear conductivity is solved exactly also in the limit of 
absence of e-e interaction (see e.g. [18] ) . Using Eq. (8) , we 
find an = n e e 2 T c bi/m 1 b x = 1/2, b y = 3/2. Thus, in this 
case <J vy /a xx — 3. Hence, for example, for this flashing 
cuts model the ratio <J vv /a xx = 3 is changed significantly 
when the temperature is changed from T < To to T > Tq. 



Quadratic response 



Here in Tj index i is axis index (x or y). For the consid- 
ered systems from the relations (6)- (8) we have r» = f c /ai 



The stationary ratchet flow appears in the second order 
of ac-force F. In this case we can operate in a similar 
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way as before. The nonlinearity occurs due to the field 
term in the kinetic equation and nonlinear e-e collision 
operator: 



d t fW-(i el + i' ee )(f {2) ) = 
-(F(t)d p )fW + i': e (f (1) *f {1) )- 



(16) 



The projection of anisotropic elastic collision opera- 
tor onto the vector functions kills the third rank ten- 
sor needed for photogalvanic current. So inclusion of 
anisotropy should be done a bit more accurately. In short 
the stationary ratchet current is generated in a follow- 
ing way. The oscillating distribution function with vec- 
tor anisotropy is converted by nonlinear e-e interactions 
to the static second angular harmonics which in turn is 
partially suppressed by linear e-e interactions and then is 
transformed to the static vector anisotropy by anisotropic 
elastic collisions. The main contribution to the station- 
ary flow reads 



Ji 



iReV^/- 1 ^)/- 1 ^)^*/^), (17) 



S 

p 



The Eq.(17) has simplified form in accordance with the 
smallness of the elastic antisymmetric operator as 
compared with the inelastic scattering obligatory 
contains higher angular harmonics). The subsequent sim- 
plifications include: the substitution of I' ee instead of I'J e , 
according Eq. (4); use of the fact that inverse operators 
J -1 do not contain antisymmetric operators; the cancel- 
lation of acting on the second angular harmon- 
ics; the replacement of the left operator (taking into 
account summation with Vi) by the inverse projected op- 
erator. As a result, we arrive at 



Vf,i = -|c , ^a J - H r i Re(r^T* fe F^F^ fe ). 



(18) 



Here e = mv 2 /2 is the particle energy, a.ijk =< 
ViVjl^Vk > t c /v 3 is a numerical tensor, characterizing 
the asymmetry of scatterers (< ... > stands for average 
over angles in the momentum space), prime again means 
the derivative over particle energy, l/r wi = — iu + l/n. 
For the specific cases of our models we obtain 



1 

48' 



(for static cuts), (19) 



l xyy = J2' ( for scmi - disks )> 



(for flashing cuts). 

07T 



For C we have: 



C 



/o°° &(/<">') V/r e )' 

J °° desifoyy ■ 



In the case of static cuts or semi-disks C = (3/2 + 
s)vp/(t c (ef)£f) (strongly degenerate Fermi case) and 



C = 24(3/2 + s)/(2 1 / 2+s mT c (T)Vrh 1 T) (Boltzmann 
case); s = 1/2, d s = 1 for fixed obstacles and s = 0, do = 
y/n/2 for flashing cuts (in this case r c (e) = const). 

To compare with results of numerical calculations it is 
convenient to write expressions for ratchet velocity com- 
ponents. For the linear polarization of monochromatic 
force we obtain 

Vf x /vT = -B(Fv T f c /T) 2 a xxx x 

[cos 2 6/(4(1 + JtD) - sin 2 9/(a x al (1 + c 2 r 2 ))] 

V fy/ V T = -B(Fv T f c /T) 2 a xyy x 

sin(20)(l + LU 2 T x T y )/(a x a 2 y (l + u 2 r 2 )(l + u 2 r 2 )) , (20) 

where B = CT f c /2v\ and we remind that Tj = r c /a,. 
For the flashing cuts model we have t c = r c , a xxx = 
1/67T, a xyy = -1 /37T, a x = 3/2 + 4/n 2 , a y = 1/2, 
C = 2^/n/(2m,T c \ / 2mT), t c = const and for the semi- 
disks model a xxx = ~a xy y — 1/12, a x = 2/3 + 8/37T, 
a y = 2/3, C = 2/(mT c Vn^T), f c cx T" 1 / 2 . Here we give 
the results for the Boltzmann distribution /(°), but sim- 
ilar calculations work for other / < -°- ) , e.g. for the Fcrmi- 
Dirac distribution. We also give a simplified derivation 
of the ratchet flow in the Appendix. It is based on the 
local equilibrium distribution and give the same results 
as Eqs. (20). 

The opposite limit in absence e-e interactions was ana- 
lytically studied for the cases of static cuts [16] (exactly) 
and approximately, for weak anisotropy, for static cuts 
or semi-disks [18]. It is important to emphasize that in 
both limits of weak and strong e-e interactions the cur- 
rent does not contain the strength of interactions. The 
transition between the regimes occurs when the interpar- 
ticle scattering rate becomes comparable with the rate 
elastic scattering on anti-dots and impurities. 



IV NUMERICAL RESULTS 

For the flashing cuts model the theory (20) gives a 
good description of numerical data (see Fig. 4). For the 
semi-disks model the agreement between the theory and 
numerical simulations (Fig. 3, top) is less accurate, e.g. 
4th ^-harmonic for Vf^ x is absent in (20). To understand 
the origins of this difference we present the map of local 
flow velocities at various polarizations 6 in Fig. 5. For 
9 = the results clearly show the appearance of turbu- 
lent flow with two vertexes behind the semi-disk. When 
the interaction scattering time tk is increased by a factor 
25 the interaction and turbulence practically disappear 
and the average local flow becomes laminar (see Fig. 5 
top left and bottom right panels). At the same time even 
with strong interactions the flow has much more laminar 
structure for 9 = n/4 (Fig. 5, top right panel) when the 
absolute value of the total ratchet velocity has its max- 
imal value (see Fig. 3 top panel). Thus the ratchet flow 
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FIG. 4: (Color online) Comparison between theory (20) (full 
curves, no adjustable parameters) and numerical data for in- 
teracting particles in the flashing cuts model (symbols) ; circles 
are for Vf yX and 9 = n/2 (here Vf tX > 0), diamonds are for 
\vj x \ and 9 = (here Vf lX < and we use the absolute value 
of Vf tX in the ratio Vf lX /vf, y (ir/4:)); Vf, y is taken at 9 = 7r/4; 
other parameters are as in Fig. 2, top panel. 
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FIG. 5: (Color online) Map of local averaged velocities in 
(x/R, y/R) plane of the semi-disks model for parameters of 
Fig. 3 (top panel) at 9 = (top left); 9 = tt/4 (top right); 9 = 
7r/2 (bottom left); 9 = and 25 times increased interaction 
time compared to other panels (tk/th = 0.5, trvt/tii = 4.7 
point in Fig. 6). The velocities are shown by arrows which 
size is proportional to the velocity amplitude, which is also 
indicated by color (from yellow/gray for large to blue/black 
for small amplitudes). 



of interacting particles has certain similarities with a hy- 
drodynamic flow of the Navier-Stocks equation around 
semi-disk body [26]. However, for 9 = ir/2 the ratchet 
flow is composed from two alternative flows moving in 
opposite directions at the cell boundaries and the semi- 
disk center (Fig. 5, bottom left panel), such a rather flow 
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FIG. 6: (color online) Dependence of the ratchet velocity v/ 
on the Kapral interaction scattering time tk in the semi-disk 
model, numerical data are shown by symbols: Vf tV /vT (red 
circles) and vj iX /vt (yellow diamonds) for 9 = n/4; vj, x /vt 
(violet squares) for 9 = 0; vj\ x /vt (blue triangles) for 9 = n/2; 
other parameters are as in Fig. 3, top panel, curves are drown 
to adapt an eye. 



is different from hydrodynamic flows with fixed velocity 
far from the body. For a qualitative description of the 
turbulent flow we may argue that the turbulence leads 
to a difference of pressures on different sides of the scat- 
terer producing different resistances for different flow di- 
rections. This generates the ratchet flow for the ac-force 
driving. In general the kinetic description is applicable 
when the interaction scattering length is large compared 
to the scatterer size, e.g. vttk > I'd for semi-disks. At 
small values of tk this condition is broken (Figs. 3,5) 
and we have transition to the hydrodynamic like regime 
where the theory (20) gives only approximate descrip- 
tion. For the flashing cuts model the kinetic description 
remains always valid since the size of scatterer is zero. 

The dependence of the ratchet velocity on the interac- 
tion scattering time tk is shown in Fig. 6. The increase 
of interactions (small tk) can change the sign of the flow 
in x-direction that is in a qualitative agreement with the 
theory (20). For weak interactions the flows are oppo- 
site in x for polarization 9 = and 9 = n/2 while at 
strong interactions they are collinear. Thus in 2DEG 
in AlGaAs/GaAs heterojunctions where interactions are 
relatively weak (r s ~ 1) the flows are opposite for two 
polarizations in agreement with the experiment [19], but 
for other materials with stronger interactions (e.g. SiGe 
with r s w 6) the flows may become collinear. We also 
note that at strong interactions the rescaled ratchet char- 
acteristics are not sensible to the temperature variation 
that indicates that we have an effective liquid flow with 
temperature independent viscosity. 
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V DISCUSSION 



Maxwell distribution ./o( v ) wc obtain 



In conclusion, our extensive numerical simulations 
show that even in the regime of strong interactions 
between particles a stationary ratchet flow is gener- 
ated by monochromatic driving in the asymmetric pe- 
riodic arrays. The obtained result are well described by 
the analytical theory based on the kinetic equation for 
strongly interacting particles. It is interesting to note 
that for asymmetric arrays the tensor of conductivity be- 
comes temperature dependent due to interplay of inter- 
actions and relaxation of high momentum harmonics (see 
Eq. (15) and discussion there). It would be interesting to 
investigate the effects of interactions on ratchet transport 
in experiments similar to those of [19]. 
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MICONANO and NANOTERRA (France) and RFBR 
NN 08-02-00152-a, 08-02-00506-a (Russian Federation). 



APPENDIX 

Here, on the example of the flashing cuts model we give 
a more simple and heuristic derivation of the ratchet flow 
compared to the exact kinetic equation approach (1). In 
the regime of very strong interactions we can assume that 
the ensemble of particles is in a local equilibrium state 
and hence the distribution function can be written as 



/(v,t) = / (v-v (t)) 



(A-l) 



where v (t) is the instantaneous velocity of the center 
of mass. We put here the particle mass m = 1. It is 
assumed that the interactions give rapid relaxation to the 
local equilibrium distribution / (v — v (i)). The matrix 
of conductivity can be determined from the momentum 
balance between acceleration created by a small applied 
static force F and momentum loss on the asymmetric cut 
scattercr: 

= F - - (< V c (v)/o(v - v ) > -v ) = , (A-2) 
at t c 

where V c (v) is the vector of average velocity after a scat- 
tering with the incident velocity v. It is expressed via the 
scattering probability W(v', v) (see the main text above) 
as V c (v) = J dv'W(v' ,v)v' . This gives the relation 



< V c (v)/o(v- v ) >= 

v y 9(v x ) 



(A-3) 



/o(v - v ) 



where 9(v) is the Heaviside function. In the linear re- 
sponse regime we can expand in Vq that gives /o(v — 



v ) = /o(v) + /o(v) : 



< V c (v)/o(v- v ) >= 

VyO(v x ) 



(A-4) 



V0,x 



^r 2 

2 v o,y 

where v is the velocity of the stationary flow. Then the 
moment balance gives 



V ,y 



■F 

1 X 



2F„ 



(A-5) 



and therefore a yy ja xx = t v /t x = a x /a y = 3 + 8/n 2 that 
is in agreement with the kinetic theory result given in the 
main text. It is interesting to note that for the nonintcr- 
acting particles we have a yy /a xx = 3 (see [18]). 

To compute the ratchet flow we should expand the lo- 
cal velocity in Eq. (A-l) up to the second order in the 
driving force F: v (t) = fjF(f) +Vf , Vf = 0(F 2 ) where 
Ti are the above values t x , T y given by the linear response; 
we note that second frequency harmonics e ±t2uJt are elim- 
inated by the time averaging. Then the time averaged 
distribution function is 



/(v) =</(v,t) > t = 

/o(v) + /o(v)-=- + 



(A-6) 



4T 2 



/o(v) 



where /o(v) = exp(— ^) is the Maxwell distribution 
and v F is the average ratchet flow velocity. Again, the 
time averaged momentum balance equation reads 



(<V c (v)/(v)>-v,)=0. 



(A-7) 



Using Eq. (A-7) we obtain from the second term the 
contribution < V c (v)/ (v) ^ >= ((1 - a x )v fx ,(l - 
a y )vf y ) which is similar to the linear response term. 
The integration of the third Gaussian term gives 
the additional contribution (F 2 /8\/2~^T)[(-t 2 cos 2 9 + 
T 2 sin 2 9) 7 2t x t v sin(2#)]. Finally we obtain 



Vfx 

v fv 



F 2 



-T? COS 2 I 



t x t 2 sin 2 i 



2t x t 2 sm{29) 



(A- 



For lo — > these expressions are in agreement with 
Eqs. (20) obtained by the kinetic equation theory. 



. After integrating over the 
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